Visualization of superposition of macroscopically distinct states 
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We propose a method of visualizing superpositions of macroscopically distinct states in many- 
body pure states. We introduce a visualization function, which is a coarse-grained quasi joint 
probability density for two or more hermitian additive operators. If a state contains superpositions 
of macroscopically distinct states, one can visualize them by plotting the visualization function for 
appropriately taken operators. We also explain how to efficiently find appropriate operators for 
a given state. As examples, we visualize four states containing superpositions of macroscopically 
distinct states: the ground state of the XY model, that of the Heisenberg antiferromagnet, a state 
in Shor's factoring algorithm, and a state in Grover's quantum search algorithm. Although the 
visualization function can take negative values, it becomes non-negative (hence becomes a coarse- 
grained joint probability density) if the characteristic width of the coarse-graining function used in 
the visualization function is sufficiently large. 

PACS numbers: 03.65.-w, 03.65.Ta, 67.40.Db, 03.67.Lx 



I. INTRODUCTION 

Visualization functions, such as the Wigner distribu- 
tion function and the Husimi function Q , are very use- 
ful. By plotting them, one can visualize quantum states 
to understand structure of these states. Furthermore, 
because these functions are, in some senses, probability 
densities, one can interpret various experimental results 
by using these functions 0. Although there are many 
methods of visualizing quantum states with small degrees 
of freedom |3|, those of visualizing quantum many-body 
states are very few 0, 0, El ■ It is therefore important 
to develop methods of visualizing quantum many-body 
states. 

In quantum many-body systems, which include quan- 
tum computers [H H 0, UJ, there are many states 
which contain superpositions of macroscopically distinct 
states 

El E M El ES E H El M M M M 

I25L l2?| . Existence of a superposition of macroscopically 
distinct states in a many-body pure state can be identi- 
fied by an index p (1 < p < 2) |2j 111 111 111 |27| : If a 
given pure state has p — 2, it contains a superposition of 
macroscopically distinct states [^J |24j . 

If every macroscopic superposition could be reduced 
to an equal-weight superposition of two macroscopically 



distinct states. 



such as |0 ■ 



■0) + 



1 • • • 1), visual- 



ization of macroscopic superpositions would be a triv- 
ial task. However, there are many other states in which 
many macroscopically distinct states are superposed with 
various weights jU III IIHH lU . Therefore it is 
also important to develop good methods of visualizing 
such complicated superpositions. 
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In this paper, we propose a method of visualizing su- 
perpositions of macroscopically distinct states contained 
in states having p = 2. We first introduce a function 
S(Ai, • ■ • , A m ), which is interpreted as a coarse-grained 
quasi joint probability density for hermitian additive op- 
erators A\, - ■■ ,A m . We next explain how to find ap- 
propriate A\, - ■ ■ ,A m efficiently for a given pure state. 
One can visualize superpositions of macroscopically dis- 
tinct states contained in a given pure state having p = 2 
by plotting 2(Ai, • • • , A m ) for appropriate A\, ■ ■ ■ , A m . 
As examples, we visualize four states having p = 2: the 
ground state of the XY model, that of the Heisenberg an- 
tiferromagnet, a state in Shor's factoring algorithm 
and a state in Grover's quantum search algorithm . Al- 
though 2 can take negative values, like the Wigner distri- 
bution function, it becomes non-negative, hence becomes 
a coarse-grained joint probability density, if the charac- 
teristic width of the coarse-graining function used in 2 is 
sufficiently large. 

This paper is organized as follows. After briefly re- 
viewing the index p in the next section, we introduce 2 
in Sec. IIII Al and explain how to find appropriate oper- 
ators efficiently in Sec. IIII Bl We visualize four states in 
Sec. IIVI Discussion is given in Sec. E] 



II. INDEX p 

To establish notation, and for the convenience of the 
reader, we briefly review the ind ex p in this section. For 
details, see Refs. [H E il |H Hi ^ • 

We first fix the energy range of interest. It determines 
the degrees of freedom of an effective theory which de- 
scribes the system under consideration. We assume that 
the system is, in that energy range, described as an A^- 
site lattice. Throughout this paper, we assume that N is 
large but finite. 
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For simplicity, we here consider only pure states, al- 
though the definition of superposition of macroscopically 
distinct states has been successfully generalized to mixed 
states |25|. Furthermore, we assume that states are 
spatially homogeneous, or effectively homogeneous as in 
quantum chaotic systems |27j or in quantum computers 
[22ll23| . For such states, we can consider a family of sim- 
ilar states {\iP N )}n- For example, each member of the 
family {\Eq)}n of the ground states of the XY model 
is the ground state \Eq) of the XY Hamiltonian of an 
TV-site system. 

The index p is defined for such families of similar states. 
For simplicity, we represent a family of states {\4> N )}n 
by a representative state (= \i/j n )). The index p 
(1 < P < 2) of \tp) is then defined by 



max(i/;|(A^) 2 |V) = 0(N P ), 

A 



(1) 



where A A = A— {i()\A\ip) |2q|. and the maximum is taken 
over all hermitian additive operators A. Here, an additive 
operator A is a sum of local operators: A = J2i=i °(0> 
where a(l) is a local operator, which is independent of 
N, on site I. We do not assume that a(l') (V ^ I) is the 
spatial translation of a(l). 

If p = 2, there is a hermitian additive operator which 
'fluctuates macroscopically' in the sense that the relative 
fluctuation does not vanish in the limit of N — ► oo: 



(Vl(Ai)^) 



N 



(N-kx). 



(2) 



Because \tp) is pure, the reason for the macroscopic fluc- 
tuation is that eigenstates of A corresponding to macro- 
scopically distinct eigenvalues are superposed with suf- 
ficiently large weights in \ip). Here, two eigenvalues 
A and A' are macroscopically distinct if and only if 
A — A'~ O(N). Therefore a pure state having p = 2 con- 
tains a su perp osition of macroscopically distinct states 
(see Refs. |24 Elf for detailed discussion). On the other 
hand, if p < 2, all additive operators 'have macroscop- 
ically definite values' in the sense that relative fluctua- 
tions of all additive operators vanish as N — > oo. In 
this case, there is no superposition of macroscopically 
distinct states in In short, one can judge whether 
a pure state contains a superposition of macroscopically 
distinct states or not by calculating the index p. 

There is an efficient method of calculating p. For sim- 
plicity, we henceforth assume that each site of the lattice 
is a spin- 1/2 system. For a given pure state \if>), we define 
the variance-covariance matrix (VCM) by 



Vai,pi> = (i>\Aa a (l)Aap(l')\i>), 



(3) 



where a, (3 = x, y, z; 1,1' = 1, 2, ■ • • , N; a x (l), & y (l), and 
<7 2 (7) are Pauli operators on site I. The VCM is a 3Nx3N 
hermitian non-negative matrix. If e\ is the maximum 
eigenvalue of the VCM, then e\ = 0(N P ~' 1 ), as shown 
in Appendix 1X1 One therefore has only to evaluate e\ to 
calculate p. 



III. VISUALIZATION METHOD 

By calculating the index p, one can judge whether 
a pure state contains a superposition of macroscopi- 
cally distinct states or not. From p only, however, one 
cannot know detailed structures of the superposition of 
macroscopically distinct states, including which macro- 
scopically distinct states are superposed and with what 
weights they are superposed. In this section, we propose 
a method of visualizing these structures of superpositions 
of macroscopically distinct states. 



A. Visualization function S 

Let A — J2i o(0 and B = J2i &(0 be hermitian additive 
operators. We assume that [A, B] ^ 0, so that the joint 
probability distribution for A and B does not exist in 
general. 

For macroscopic systems, one is usually interested in 
states in which typical values of additive operators are 
0(N) 28]. Typical values ofA/N and B/N are therefore 
O(N ). Their commutator is small in the sense that 



A B 

N ' N 



(4) 



because [a(l),b(l')] = for I ^ I' [H HI H|. Since 
N is large but finite, the above commutator does not 
vanish. In real experiments, however, resolutions of mea- 
surements are limited. Equation Q therefore indicates 
that noncommutativity of additive operators could not 
be detected for large V. This suggests that we may be 
able to introduce a function which can be well regarded 
as a coarse-grained joint probability density for A and 
B. Note that the finite resolution is essential, because 
noncommutativity, however small it is, can be detected 
if the resolutions of experiments are high enough [2{j . 

We formulate the above idea as follows. Consider the 
spectral decompositions of A and B: 



A = E av a( a i b=y, bv b( b ) 



(5) 



age i 



BeE f 



where and are the spectra of A and B, respec- 
tively, and Vt(A) and V&{B) are the projection oper- 
ators onto the eigenspaces of eigenvalues A and B, re- 
spectively. To take account of finite resolutions of exper- 
iments, we smear the projection operators to obtain 



r A (A)^ Yl ^a(a,a')p a (a% 



(6) 



A' EE > 



and similarly for Vg(B). Here, A is a real continuous 
variable (A S R), and ws(A, A') is a coarse-graining 
function. It centers at A = A' with a characteristic width 
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Wx, and satisfies 



wt{A,A')>Q for all A, A', 



(7) 



/+oo 
w A (A, A')dA = 1 for all A'. (8) 
-oo 

The coarse-graining functions w A ,Wg should not have 
complicated forms; they should be physically reasonable 
ones. To be definite, we henceforth assume that Wx = 
W & = W and w A {X,X') = w § (X,X') = w(X - X'), 
where 



w(X) = 



1 



2nW 



exp 



X 2 \ 



2W 2 



(9) 



Clearly V A {A) and "P g(B) are non- negative hermitian 
operators satisfying 

/+oo r + oo 

V A {A)dA = / V § (B)dB = 1. (10) 
-oo J —oo 

They give coarse-grained probability densities 2^ and E ^ 
for A and B, respectively, for a given pure state 1^) by 

E A (A) = MV A (A)\i>) (4 6 1), (11) 
R A {B) = (ij\V e (B)\^) (B e R). (12) 
Now we define 

E(A,B) = ^\V A (A)V 6 (B)+V^B)r A (AM (13) 

for A, B E R. One can easily verify the following: 



(A, B) is real, 

E(A, B)dAdB = 1, 



(14) 
(15) 



/-t-oo /•+OC 
E(A, B)dB = E A (A), / E(A,B)dA = E$(B).(16) 
-oo J —oo 

In general, E(A, B) can take negative values. If it is 
non-negative, Eqs. Ijl5|l and l|16(l show that it can be 
interpreted as a coarse-grained joint probability density 
(cgJPD) for A and B. In fact, as we will demonstrate in 
the following sections, E(A, B) becomes non-negative if 
W and N are large enough, for many states of interest. 
Furthermore, even if W and N are not large, negative- 
valued regions of E(A, B) are small. In this case, E(A, B) 
can be considered as a coarse-grained quasi joint proba- 
bility density (cgQJPD) for A and B. 

The non-negativity of E(A, B) becomes obvious as 
W -> oo, for which E(A, B) ~ w(A)w(B) > for all 
A, B. For smaller W, the smallest value of W that makes 
E(A, B) non-negative depends on A, B and \ip). There- 
fore, in general, the non-negativity should be checked a 
posteriori. 

We can also introduce E for m (> 3) hermitian additive 
operators A% , A2, ■ ■ ■ , A m by 

E(A U --- ,A m ) 

= i E^l^ (1) (4r(l)) ' ' -^ (m) (iw)l*. (17) 



where the sum is taken over all permutations n of the 
numbers 1, 2, • • • , m. 

If 1^) has p — 2, one can visualize structure of the 
macroscopic superpositions contained in by plotting 
E(Axy ,A m ) versus (At,--- ,A m ), if ii,--- ,i m are 
appropriately taken. We call such a plot a visualization of 
superpositions of macroscopically distinct states in 
An efficient method of finding appropriate A\, - ■ ■ ,A m 
will be explained in the next subsection. 



B. Efficient method of finding appropriate 
operators 

In principle, one can take any hermitian additive op- 
erators Ai,--- ,i m , and plot E(Ai,--- ,A m ). In this 
paper, however, we are interested in states having p = 
2, which contain superpositions of macroscopically dis- 
tinct states. Such superpositions are characterized by 
macroscopic fluctuations of certain additive operators 
(see Sec. Inland Refs. [U Q). Therefore, as will be 
demonstrated in the next section, we can visualize such 
superpositions by including macroscopically fluctuating 
operator(s) in ii,--- , i m of E(Ax,--- ,A m ). In this 
subsection, we present an efficient method of finding a 
set S of macroscopically fluctuating hermitian additive 
operators. 

For a given pure state we diagonalize the VCM to 
obtain its eigenvalues, e\ > e2 > • • • > e^N, and eigenvec- 
tors. From the eigenvectors, we construct a complete or- 
thogonal system: {{c^}, {c 2 t }, • • • , {c£ N }} (a = x, y, z; 
I = 1,2,- •• ,N). Here, {cj z } £ C 3N is an eigenvector 
of the VCM corresponding to . We assume that each 
cj; is asymptotically independent of N, and that we can 
normalize {cj;} as J2 al \c^i\ 2 — N. By taking an ap- 
propriate limit of {cj;} as described in Appendix iBl we 
obtain a vector {c l al }, whose elements are independent of 
N . From this vector, we construct the additive operator: 



N 



(18) 



1=1 a=x,y,z 



As shown in Appendix [U] Ai fluctuates macroscopically 
if and only if a = O(N). 

If ei = O(N) and Ai is hermitian, we let Ai be an ele- 
ment of S. If ei = O(N) and Ai is non-hermitian, on the 
other hand, we decompose A t into the real and the imag- 
inary parts: i t = if + iif 1 , where if = (i; + ij)/2 
and i™ = (Ai~Al)/2i. It is known that if and/or i™ 
fluctuate(s) macroscopically (see Ref. |24( and Appendix 
A). We let such macroscopically fluctuating part(s) be an 
element(s) of S. In this way, we obtain a set of macro- 
scopically fluctuating hermitian additive operators, e.g., 
as S = {A x , A^ , A z , A T £ , A\ m , As) . Several examples of 
S will be given in the next section. Any macroscopi- 
cally fluctuating additive operator includes at least one 
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element of S as a component in the sense explained in 
Appendix ID1 

The number of the elements of iS is O(N ), because 

d > and J2 3 i=i ei = Eili T, a =x, v ,z V <xl,<xl < 3N. One 
can obtain S efficiently, because one has only to diago- 
nalize the VCM, which is a 3iV x 3N hermitian matrix. 

By including an element(s) of S into A\,--- ,A m 
of E(A\,--- ,A m ), one can visualize superpositions 
of macroscopically distinct states in \ip) by plotting 
E(A U --- ,A m ). 

IV. EXAMPLES 

To demonstrate usefulness of the visualization method, 
we visualize four states having p = 2 in this section. 

A. XY model 

First, we visualize the exact ground state of the XY 
model on a two-dimensional square lattice of N sites. 
The Hamiltonian is 

h = -J2 ^*(o^(0+^(o^(0], as) 

<M'> 

where < 1,1' > denotes the nearest neighbors. If N is 
finite, 'ground states' obtained by the mean- field ap- 
proximation are very different from the exact ground 
state jaSEJESEl- These mean-field ground states 
are degenerate symmetry-breaking states with non-zero 
order parameters. They are separable states, because the 
mean-field approximation neglects the correlations be- 
tween sites. On the other hand, the exa ct g round state 
is unique, symmetric, and has p = 2 p3. 130L l3ll l32l l33l| . 
We visualize the exact ground state. 

By numerical calculations, we find that e\ — e 2 = 
<D(N), e t = o{N) (i > 3), ii = EjliMO = M x , and 
A 2 = Ei=i Oy(0 = My. Hence, 

S = {M x ,M y }. (20) 

In Fig. □ we plot E(M x ,M y ) for N = 14 without 
coarse-graining, i.e., W ^ 0, for which the coarse- 
graining function of Eq. becomes the delta function 
S(X). In the figure, cS(0) [c € K] is represented by a ver- 
tical line with height c. Positive values are represented 
by red vertical lines, whereas negative values are repre- 
sented by blue vertical lines. Because S(M x ,M y ) takes 
negative values at some points, it is not a JPD for M x 
and M, y . 

However, negative values are expected to approach 
as W is increased. To see this, we plot in Fig. [2] the 
integral 1^ ^ of negative values versus W, where 1^ g 
is defined by 

j /;;»»™. (21) 



It is seen that 1^ ^ indeed approaches as W is in- 
creased. H(M X! My) therefore becomes a cgJPD if W is 
sufficiently large. 

For example, we plot E(M x ,M y ) with W = 3.2 in 
Fig. |3 In this case, E.(M Xl M y ) is non-negative, and 
therefore it is a cgJPD. From this figure, we can clearly 
understand the structure of the superposition of macro- 
scopically distinct states: Many macroscopically distinct 
states which have macroscopically definite U(l) order pa- 
rameters (M x ,M y ) are so superposed that the ground 
state has the U(l) symmetry. 

In Fig.H we plot E(M x ,M y ) with smaller W, W = 2. 
Because E(M X , M y ) has negative- valued regions, it is not 
a JPD. However, since A ~j | is small, E(M x ,M y ) is 
well regarded as a quasi JPD. The figure shows the same 
U(l)-symmetrical structure as that of Fig- EI 

One can utilize either Fig. |31 or Fig. 0] depending on 
the purpose: When one wants a cgJPD, Fig.EJshould be 
used. On the other hand, when one wants to see more 
detailed structures, including quantum effects that make 
S negative, then Fig. 0] (or, Fig.Q would be better. In 
this way, one can adjust W to obtain a useful 5 according 
to the purpose. 



XY, N=14, W-->0 
negative 




FIG. 1: (Color online) E(M X ,M V ) with W -> for the exact 
ground state of the XY model on a two-dimensional square 
lattice with TV = 14. c5(0) [c £ R] is represented by a vertical 
line with height c. Positive values are represented by red 
vertical lines, whereas negative values are represented by blue 
vertical lines. 



B. Heisenberg antiferromagnet 

Second, we visualize the exact ground state of the 
Heisenberg antiferromagnet on a two-dimensional square 
lattice of N sites. The Hamiltonian is 

H= [**(0** (O + *v(0*v(O+ *«(0*»(O] -(22) 
<l,l'> 
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FIG. 2: The integral A - f of negative values of S(M X , M y ) 
versus VF for the exact ground state of the XY model on a 
two-dimensional square lattice with N = 14. 



FIG. 4: (Color online) E(M x ,M y ) with W = 2 for the exact 
ground state of the XY model on a two-dimensional square 
lattice with N = 14. Positive-valued regions are colored red, 
whereas negative- valued regions are colored blue. 



XY, N=14, W=3.2 

0.0016 



S(M X ,M ) 
0.002 
0.001 





FIG. 3: (Color online) a(M x , M y ) with W = 3.2 for the exact 
ground state of the XY model on a two-dimensional square 
lattice with N = 14. 



The 'ground states' obtained by the mean-field approxi- 
mation are degenerate, symmetry-breaking, and separa- 
ble. On the other hand, the exact ground state is unique, 
symmetric, and has p = 2 if N is finite |2ll3lLl3ll3ll3^. 
We visualize the exact ground state. 

By numerical calculations, we find that e\ = 62 = 
O(N), e, = o(N) (i > 4), Ai = E^i(-l)'5-*(0 = 
M = Eili(-1)'^(0 = K, and A 3 = 
^-ly&zQ) = Mf. Hence 



63 = 



14. Because 3(M** , My) is non-negative, it is a JPD for 

MJ* and M*'. It is also seen that many macroscopically 
distinct states are superposed in the ground state. 

By increasing W , we obtain more understandable pic- 
tures. For example, in Fig. we plot 3(M£ t ,MJ t ) with 
W = 2. It is seen that many macroscopically distinct 
states are so superposed that the ground state is sym- 
metric, like the ground state of the XY model. 



HAF, N=14, W->0 



E(M x st ,M, * 




FIG. 5: (Color online) E(Mf,Mf) with W -> for the 
exact ground state of the Heisenberg antiferromagnet on a 
two-dimensional square lattice with N = 14. e<5(0) [c G R] is 
represented by a vertical line with height c. 



S = {M?,M?,M?}. 



(23) 



E(Mf,M^,Mf) is hard to plot, we plot 
Because of the rotational symmetry of the 



M st 

v 



Because 

model, E(M^,M^) = E(M^,Mf) = S(Mf,Mf ). 



In Fig. we plot 3(M* t , M**) with J4 7 — > for TV 



C. Shor's factoring algorithm 

Third, we visualize a state in Shor's factoring algo- 
rithm |3 EI ■ Let / be an integer to be factored. We use 
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HAF, N=14, W=2 




FIG. 6: (Color online) 3(Mf ,Mf ) with W = 2 for the ex- 
act ground state of the Heisenberg antiferromagnet on a two- 
dimensional square lattice with N = 14. 



two quantum registers, the first and the second registers, 
which are composed of Ni (2 log 2 I < N\ < 2 log 2 7+1) 
and iV 2 (log 2 7 < N2 < log 2 I + 1) qubits, respectively. 
We denote the total number of qubits by N = Ni + iV 2 . 
If the order r is 6, for example, the state 



and My . There are four peaks, which represent a super- 
position of approximately four macroscopically distinct 
states. 

In Fig.Hni we plot E(M x ,Mf) with smaller W,W=1. 
In this case, there is a negative-valued region. However, 
because |7^ ^ st \ is small, E(M X , My 1 ) is interpreted as 

a cgQJPD. E(M x ,My t ) again represents four peaks. We 
have also observed such a four-peak structure for some 
other values of (7,x)'s. 



Shor, (21,5), W->0 




/2JVT 



J2 K mod/) 2 |a)i, 



(24) 



a=0 



which appears just after the modular exponentiation, has 
p = 2 I22u23j|. Here, | • • • )i an d I " ")2 represent the first 
and the second registers, respectively, and x (x < 7) is a 
randomly taken integer coprime to 7. 

For the states of r — 6, we numerically find that e% — 

e 2 = O(N), e t = o(N) (i > 3), A 1 = y^f <7 X (Z) = 

M«, and A 2 = yfiY^M- 1 ) 1 ^ 1 ) = ( see 
Appendix B). Here, the qubit states are |0) and |1) 
(a z |0) = -|0), a z \l) = |1)). Hence 



S = {MU ) MfV}. 



(25) 



Because M x ^ and My H± ' fluctuate macroscopically, M x 
and My 1 also fluctuate macroscopically [23| • We here use 
M x and instead of AfP and My t{1) . 

In Fig. □ we plot E(M x ,Mf) with FT -> for 
(7,x) = (21,5). Because ^(M^M**) takes negative val- 
ues at some points, it is not a JPD. To see the behavior 
of negative values, we plot in Fig. El the integral 7,> ,- /3t 

X 1 y 

versus W. We see again that 7vv A ~, 3t approaches as W 

x ' y 

is increased. H(M S ,M*') therefore becomes a eg JPD if 
W is sufficiently large. 

In Fig. El we plot E(M x ,Mf) with W = 1.4. Be- 
cause E(M X , M^) is non-negative, it is a cgJPD for 



r st(l 



FIG. 7: (Color online) HfM,, Mf ) with W -> for the state 
just after the modular exponentiation with (I,x) = (21,5). 
c5(0) [c 6 R, c > 0] is represented by a vertical line with 
height c. c<5(0) [c g R, c < 0] is represented by a vertical 
line with height 10c, in order to make negative values more 
visible. Positive values are represented by red vertical lines, 
whereas negative values are represented by blue vertical lines. 
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FIG. 8: The integral 7,> ^ Ist of negative values of 

X ! y 

S(M X , My) versus W for the state just after the modular 
exponentiation with (I,x) = (21,5). 
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Shor, (21,5), W=1.4 
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FIG. 9: (Color online) E(M X , Mf) with W = 1.4 for the state 
just after the modular exponentiation with (I,x) = (21,5). 



Shor, (21,5), W=1 
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FIG. 10: (Color online) S(M m ,Mf) with W = 1 for the state 
just after the modular exponentiation with (I,x) = (21,5). 
Positive-valued regions are colored red, whereas negative- 
valued regions are colored blue. In order to make negative 
regions more visible, negative values are multiplied by 10. 



D. Grover's quantum search algorithm 

Finally, we visualize a state in Grover's quantum 
search algorithm |j| Let us consider the problem 

of finding a solution to the equation f(x) = l among 2 N 
possibilities, where f(x) is a function, / : {0, 1, • • • , 2 N — 
1} — * {0,1}. These 2^ possibilities are indexed by 2 N 
computational basis states, which are tensor products of 
|0) or |1) of Aqubits. Here, <7 Z |0) = -|0) and<7 2 |l) = |1). 
Let \Gk) be the state which appears after k Grover iter- 
ations. It was shown that if the number of the solutions 
is O(N ), \G k )'s whose k satisfies 



S < 



Ak + 2 



< 7T 



(26) 



have p = 2, irrespective of which numbers are the solu- 
tions |23|. Here, <5 is an arbitrary small positive constant 
being independent of N . 

To be definite, we assume that the state 11®^) indexes 
the solution. Then \Gt) is written as 



\G k 



2fc + 1 



1 



V2 W -1 



+ 



+ sin 



2k + 1 



II 



(27) 



where cos § = y^W^Tf/W, and [\0® N ) + ■ ■ ■] is the 
equal-weight superposition of all computational basis 
states except for |1® ). Among many A:'s which satisfy 
Eq. we use k = R/2 for even R, and k = R/2 + 0.5 

for odd R, where R = CI (j^ 1 arccos\/2 _Jv ^ is the num- 
ber of total Grover iterations. Here, CI(x) denotes the 
integer closest to the real number x. 

We numerically find that e\ — 0(N), = o(N) (i > 

2), and A\ = J^iLi ^^(0 + ^^(0 = ^-z (see 
Appendix B). Hence 



S = {M x . z }. 



(28) 



Because S has only one element, the macroscopic super- 
position can be visualized by plotting the probability den- 
sity (GklVgj (M x . z )\Gk)- However, because it is more 

interesting to plot Z(M x _ Zl A), where A is a hermitian ad- 
ditive operator, we plot S(M x - z ,M y ) in this paper. The 
shape of (GklV^ (M x . z )\Gk) can be deduced from that 
of 3(M K , My). °"* 

In Fig. [IB we plot E(M x _ z ,M y ) with W -> for N = 
12. Because 5 takes negative values at some points, it is 
not a JPD. To see the behavior of the negative values, we 



.M, 



versus W. I 



plot in Fig. ^] the integral I 
approaches as W is increased. E(M X . Z , M y ) therefore 
becomes a cgJPD for M x . z and M y if W is sufficiently 
large. 

In Fig.HSl we plot E(M x . z ,M y ) with W — 2. Because 
there are small negative-valued regions, it is a cgQJPD. 
It is seen that the state is approximately a cat state, 
i.e., an equal- weight superposition of two macroscopically 
distinct states. Although this information can also be 
obtained by plotting (Gk\Vj^ (M x . z )\Gk), we can see 
interesting structures of \Gk), including negative- valued 
regions, by plotting E(M X . Z , M y ). 



V. DISCUSSION 

A. Non-negativity of S 

In the previous section, we have observed that an ap- 
propriate value of W which makes 3 non-negative largely 
depends on the quantum state to be visualized. For 
example, E(M X , M y ) for the ground state of the XY 
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Graver, N=12, W->0 

negative Graver, N= 12, W=2 




FIG. 11: (Color online) E(M X . Z ,M V ) with W -> for a state 
in Grover's quantum search algorithm with N = 12. cS(0) 
[c 6 R] is represented by a vertical line with height c. Positive 
values are represented by red vertical lines, whereas negative 
values are represented by blue vertical lines. 



FIG. 13: (Color online) E(M x . z ,M y ) with W = 2 for a state 
in Grover's quantum search algorithm with N = 12. Positive- 
valued regions are colored red, whereas negative-valued re- 
gions are colored blue. 
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FIG. 12: The integral I AI $ of negative values of 
E(M X . Z , My) versus W for a state in Grover's quantum search 
algorithm with N = 12. 



model becomes non- negative with W = 3.2, whereas 
S(M X ,M**) for the state in Shor's factoring algorithm 
becomes non- negative with smaller W, W = 1.4, for 
the same value of N. Furthermore, S(MJ*,M**) for the 
ground state of the Heisenberg antiferromagnet is non- 
negative with any W. Therefore, in general, one must 
find an appropriate value of W a posteriori. 

However, it is worth mentioning that a sufficient mag- 
nitude of W which makes 3 non-negative seems to be 
O(N). To see this, consider the following three exam- 
ples. 

Example 1: In Figs. 1141 and 1151 we plot 1^ AI 
versus N for \G k ) of Eq. J23 with W = 0(N) 'and 
W = 0(y/~N), respectively. Here, k = R/2 for even i?, 



and k = R/2 + 0.5 for odd R. It is seen that 1^ A[ 
approaches as N is increased if W — O(N), whereas it 
does not approach if W = 0(V~N)- 

Example 2: In Fig. El we plot 3{M X , My) for the sep- 
arable state |0 8JV ) with W = 1.5 and N = 14. Here, 
a z \0) = — 10). There are negative- valued regions. In 
Figs. El and El we plot 1^ $ versus N with W = 

O(N) and W = 0(y/~N), respectively. We can see again 
that I A j approaches as N is increased if W = O (N) , 

whereas it does not approach if W = Q{-jN). 



Example 3: 



For the cat state 4= 

V2 



10®*) 



V2 



in which M z fluctuates macroscopically, E(M Zl M x ) and 
S(M z ,M y ) are non-negative. On the other hand, 
E(M x ,M y ) can take negative values. In Fig. El we plot 
E(M X , My) with W = 1.5 and N = 14. It is seen that 
there are negative- valued regions. In Figs. [5D| and |^ 
we plot Ife^ with W = O(N) and W = 0{VN), 
respectively. 1^ A y again approaches as N is in- 
creased if W = O(N), whereas it does not approach 
i£W = O(VN). 

From these (and some other) examples, it is expected 
that O(N) is a sufficient magnitude of W which makes 3 
non-negative for sufficiently large N. This expectation is 
reasonable, because W = O(N) means that the relative 
error of a measurement is independent of the system size 
N, which is a usual situation for macroscopic systems. 

Whether S(Ai,--- ,A m ) is non- negative or not de- 
pends also on which additive operators A\, ■ ■ ■ , A m are 
used. For the ground state of the XY model, for example, 
if we use M x and M z instead of M x and M y , S(M X , M z ) 
is non-negative with any W, because the ground state 
is an eigenstate of M z corresponding to the eigenvalue 
M z = 0, hence ~{M X ,M Z ) = ~ a1 JM x )w(M z ) > 0. 

At the time of writing, however, we do not know a 
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method of finding hermitian additive operators and W 
which make S non-negative for a given state. To find 
such a method will be a subject of future studies. 



B. Negative- valued regions of S 

If [A, B] = 0, 5(j4, B) is non-negative. It is there- 
fore expected that negative- valued regions of S represent 
some quantum natures, like those of the Wigner distri- 
bution function. 

It seems that superposition of macroscopically distinct 
states studied in this paper is not directly related to 
negative-valued regions. For example, S(M^*,M*') is 
non-negative with any W for the ground state of the 
Heisenberg antiferromagnet, which has p = 2. 

In the previous subsection, on the other hand, we have 
seen that "E,(M X , M y ) has negative- valued regions for the 
separable state |0®"). Because the separable state has 
no quantum nature other than the quantum coherence 
within each site, the negative-valued regions should rep- 
resent this quantum coherence. This expectation is rea- 
sonable, because E(M x ,M y ) is non-negative with any 
W for the random state p r = 577T, which has neither 
entanglement nor quantum coherence. Here, we pro- 
visionally define S for a mixed state p by H(j4, B) = 
\Tvp[P A {A)V h {B)+V § {B)V A {A)]. 

Detailed analysis of negative-valued regions is, how- 
ever, beyond the scope of the present paper. It will also 
be a subject of future studies. 
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The integral J A - f ^ of negative values of 
versus N with W = N/6 and W = N/8 for 
states in Grover's quantum search algorithm. Lines are guides 
to the eye. 
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FIG. 15: The integral 7. 



of negative values of 



S(M X . Z ,M V ) versus N with W = \ /f N and W = VoJN for 
states in Grover's quantum search algorithm. 
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FIG. 16: (Color online) E(M x ,M y ) with W = 1.5 for the 
separable state 10®^} with = 14. Positive-valued regions 
are colored red, whereas negative-valued regions are colored 
blue. 



APPENDIX A: ei = ©(iV^- 1 ) 

In this appendix, we show that e\ = 0(N P ~ 1 ). For a 
given pure state let {c^} S C 3N be an eigenvector of 
the VCM corresponding to the maximum eigenvalue e\. 
We normalize it as ^Z a i l^o/l 2 = ^ ■ From the eigenvector, 
we construct the operator: 



Ai 



(Al) 
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If it is hermitian and all c^'s are independent of 
N, it gives the maximum of Eq. Q). Therefore, 
max y |(?A|(AA) 2 |V>) = e±N in Eq. (JTJ. Hence e\ = 
OiNP- 1 ). 
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FIG. 17: The integral 1^ ^ of negative values of 
E(M x ,M y ) versus TV with W '= TV/6, W = TV/8, and 
W = TV/10 for separable states |0® w ). Lines are guides to 
the eye. 



FIG. 19: (Color online) E(M x ,M y ) with W = 1.5 for the 
cat state ^10®^} + -±= \1® N ) with N = 14. Positive- valued 
regions are colored red, whereas negative-valued regions are 
colored blue. 
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FIG. 18: The integral 1^ ^ of negative values of 

E(M x ,M y ) versus N with W = \/TV and W = VoIN for 
separable states 10®^}. 



If Ai is non-hermitian and all c i's are independent 
of TV, we decompose it as: A x = + iA™, where 
Af = (Ai + k\)/2 and If* = (A x - A|)/2i. Then 
(^|(Aif) 2 |V) = O(eiTV) or (^(Ai™) 2 ^) = O(eiTV), 
because 
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FIG. 20: The integral 1^ £ { of negative values of 
E(M x ,M y ) versus TV with W '= TV/6, W = TV/8, and 
W = TV/10 for cat states -±- \0® N ) + -±= H®*). Lines are 
guides to the eye. 



operator 



jAAil^H = \\AA?\i/>) + iAAT 

< HAiri^ll + IIAiri^ll. (A2) 



al 



(A3) 



Assume that (<0|(AAf ) 2 \i/j) = O(eiTV). Because Ai is 
additive, A™ is also additive. Then max x(ip\( A A) 2 \ip) = 
O(eiTV) in Eq. (QJ. Hence ei = C^TV?- 1 ). 

If some c ;'s depend on TV, we compose the additive 



where {c 1 ^} is obtained by taking an appropriate limit of 
{cjj} as described in Appendix B. It can be shown that 
(V>|AijAii|V>) = O(eiTV). In fact, by defining = 
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FIG. 21: The integral 1^ £ { of negative values of 
E(M x ,M y ) versus AT with W = and W = \/05iV for 



cat states -^\0® N ) + \1® N ) 



r 1 

c al 



-al' 



|A4AAi|V> 



„1* 1 T ,T 



= ex AT + e\o{N) 
= 0( ei N), 



aipi' 



where we have used the facts that {c ,} is an eigenvector 
of the VCM corresponding to ex, that Ic^l 2 = AT, 



that d^i 



o(N°), and that < J2 a w> d a \*d^,V aW v < 
e 1 N. If A 1 is hermitian, ma,x A (ip\(AA) 2 \ip) = 0{e 1 N) in 
Eq. Hence ei = 0(N P ^ 1 ). If Ai is not hermitian, its 
real or imaginary part gives 0(e\N) fluctuation. There- 
fore, max^KAi) 2 !^) = Q{e 1 N) in Eq. JTJ. Hence 
ex = ©(A" 5 - 1 ). 

In conclusion, we have shown that d = ©(A"'- 1 ). 



APPENDIX B: COMPOSITION OF {t&j} FROM 



of elements among 3Af elements of {c^, ; } can be modified, 
because it does not alter the leading term (with respect 
to AT) of (-0 1 A A j A ^4.^ | -0) . Using this property, we can 
adjust Ai for our convenience. 

For the state of Eq. (|24[) with r = 6, for example, 



= J y/N/jNi - 1) (a = x; 2<l< Aq), 
(others), 



-ai 



= / (-ll'^jpri) (a = y; 2 < J < Aq), 
(others). 



We therefore obtain 



/3/2 (a = z; 1 < i < Aq) 
(others), 

/3/2 (a = y; 1 < I < ATx) 
(others). 



(B2) 
(B3) 

(B4) 
(B5) 



Or, we can modify the / = 1 terms of these results as 
= in accordance with the I — 1 terms of 



al 



a 1 



Eqs. (|B2(1 and (|B3|) . We have employed the latter forms 
in Sec. UVTI 

Moreover, for \Gk) with k = R/2 (even i?) or k — 
R/2 + 0.5 (odd R), 



r 1 



-a/Va 2 + b 2 + c 2 (a = x; 1<1<N), 



ib/Va 2 + b 2 + 



b 2 



- 2 (a = y; 1 < I < N), (B6) 
? (a = z; 1 < Z < AT). 



Here, a, 6 and c are real numbers which depend on N. 
It is numerically shown that limjv_ +00 (a — c) = and 
(A4) limjv— &/V a 2 + b 2 + c 2 = 0. We therefore obtain 



-1/V2 (a = x; 1 <l < AT), 
(a = y; 1 < I < AT), 
l/\/2 (a = z; 1<1<N), 



(B7) 



which has been used in Sec. IIV Dl 



APPENDIX C: Ai FLUCTUATES 
MACROSCOPICALLY IF AND ONLY IF a = O(N) 

In this appendix, we show that Ai of Eq. I|18|) fluctuates 
macroscopically if and only if e, = O(N). Using c l al of 



Appendix El we define 



Then 



By diagonalizing the VCM, one obtains {c Q \} £ C 3N 
corresponding to ej. Each element cj ; generally depends 
on Af, whereas c^ ; 's composing the additive operator Ai 
through Eq. H18[) should be independent of AT. We can 
deduce c l al from simply as follows. 

Let us define a parameter v = l/N, and denote c Q l ; by 
c* v (N). We take the following limit: 



N 



lim c* v (N') 



(Bl) 



where ^ is kept constant in this limit. Then, c l al is given 
by & al = & a (l/N). Note that a small number [= 0(N )} 



\aA\aa\iP) 



^ C al C pl< V otl,0V 

Q//3/' 



aipi' 



pi' 



eiiV + e i0 (AT) + o(iV 2 



(CI) 



where we have used the facts that {cj z } is an eigenvector 
of the VCM corresponding to a, that J2 a i \^ai\ 2 = ^> 



and that <i. 



o(JV°). Therefore, if = O(AT) then 
|AAjA^|0> = 0(N 2 ). On the other hand, if a = 

it A A.\J,\ — „( AT2\ 



o(N) then (i/j\AA]AAi\ip) = o(N 2 ). 
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APPENDIX D: ANY MACROSCOPICALLY 
FLUCTUATING ADDITIVE OPERATOR 
INCLUDES AN ELEMENT OF S 

For an additive operator A — 5Z Q ; c Q ;0" Q (Z), the coef- 
ficient vector {c a i} <E C 3N can be expressed as a linear 
combination of {cjj's: c a i = X^&cA, where &'s are 
coefficients satisfying X^l^l 2 = O(N ). Assume that 
& = (N°) if ei = O(N) (i = 1, • • • , 3JV). Then 

(^|AitAi^) = iV^ |^ 2 e» = o(7V 2 ), (Dl) 

i 

where we have used the facts that {c Q \}'s are orthogo- 
nal eigenvectors of the VCM, and that J2 a i \^ai\ 2 = N. 
Equation (|D 1|) shows that A does not fluctuate macro- 
scopically. 



In other words, if A fluctuates macroscopically, its co- 
efficient vector {c a i} includes at least one {c s ; } whose 
ej = O(N) as a component of the linear combination 
with the weight & = G(N°). Therefore, 

i = ^[•••+&c>a(0 + --] 

= + + (D2) 

which shows that A includes Ai (hence also A™ and A 1 ™) 
with the weight £j = O(N ). In this sense, at least one 
element of 5 is 'included' in A, if j4 fluctuates macro- 
scopically. 
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